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In this work we discuss a short range version of the p-spin model. The 
model is provided with a parameter that allows to control the crossover with 
the mean field behaviour. We detect a discrepancy between the perturbative 
approach and numerical simulation. We attribute it to non-perturbative efltects 
due to the finite probability that each particular realization of the disorder 
allows for the formation of regions where the system is less frustrated and 
locally freezes at a higher temperature. 



PACS Numbers 05-7510N 



I. INTRODUCTION 

The mean field approximation is certainly a very useful starting point to study systems 
with long range interactions. Many physical features may differ from mean field behaviour 
if a model with short range interactions in finite dimension is considered. Generally speak- 
ing, the analysis of the Gaussian fluctuations is the first step to understand if the mean 
field approximation is well suited to describe the physics of the model. The perturbative 
fluctuations around the minima of the effective action may change the critical behaviour of 
the model or even destroy the phase transition. 

In spin glasses and other similar glassy models the situation is rather complicated Q. 
Generally speaking, one can distinguish two main classes of models: a first class, which com- 
prehends models that in mean field undergo a continuous phase transition with a full replica 
symmetry breaking, and a second class, comprehending models which are well described, in 
mean field, by a IRSB low-temperature solution. For this latter class, if the replica sym- 
metry is broken with a parameter qi that is of order 0(1), the transition is discontinuous. 
Models of the first class describe well the physics of real spin glasses which undergo a second 
order phase transition with a divergent non-linear susceptibility and long range correlations. 
For these models the analysis of the infrared divergences of the propagators is rather subtle 
because not all of the modes in replica space are simply related to fluctuations of a physical 
quantity |,|. 

In this work we are studying models of the second class where the transition does not 
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imply any zero mass mode. The transition is discontinuous in the sense that, for T < Tc, 
another solution, different from the high temperature one, becomes dominant in the partition 
function for large N. In this kind of transition there are no precursor effects like quantities 
that diverge as T approaches Tc from above. Physically speaking, the discontinuity of the 
transition reflects the fact that the number of the metastable states becomes non-extensive 
for low enough temperatures and so a finite probability that two replicas belong to the same 
state arises. The probability P{q) that two replicas a and b have an overlap qat — q develops 
a delta function at qi when T <Tc- The weight of this delta function goes to zero as T goes 
to Tc and it increases continuously as the temperature decreases below Tc- This implies that 
thermodynamic quantities such as the entropy are continuous at Tc and no latent heat is 
involved in the transition which therefore appears to be half way between first and second 
order. One may wonder if the strange features of this transition survive when models of the 
second class are considered in finite dimensions. 

In this work we shall introduce a short range version of the p-spin glass model. The 
model, defined in d spatial dimensions, contains an appropriate parameter M ensuring the 
mean field solution to be exact for M — > oo. 

The analysis of the perturbative fiuctuations (which are of order 0{1/M)) around the 
MF solution and of the Gaussian propagators provides some insight on the nature of the 
transition in finite dimension. A second order phase transition would imply the existence of 
zero mass modes at the critical temperature which reflect the divergence of the correlation 
length. We will find no divergence of the zero momentum propagators and therefore this 
transition seems to invoke no diverging correlations. Moreover, calculating the 0{1/M) 
corrections of the relevant thermodynamic quantities around the saddle point, we observe 
that a perturbative approach shows no qualitative change of the phase transition with respect 
to mean field: the transition still appears discontinuous and still involves no latent heat. 

Numerical simulations on the model |^-^ nevertheless indicate the existence, in d = 
3 and in d 4, of a divergent susceptibility and of all the phenomenology typical of a 
continuous second order transition at a temperature greater or equal than the mean field 
static critical temperature. 

This apparent conflict between the two results may possibly be explained considering 
non-perturbative effects: the existence of regions in space where the system is locally less 
frustrated than on average and in which the spins are long range-correlated even for some 
T>Tc. 

The paper is organized as follows. In the next section we shall briefly review the replica 
solution of the long-range p-spin glass model in absence of external magnetic field. We will 
refer to the original references for a more detailed discussion of the model and of the replica 



method. In section III we will introduce our short-range version of the model. In section 
|iv| we will calculate the Gaussian fiuctuations around the large- Af solution of the model in 
the hot and in the cold phase. In section ^ we will calculate the propagators which show no 
zero-mass modes. We will also discuss the propagators on the dynamical line i.e. obtained 
by imposing the condition of marginality. Finally, we shall draw a possible interpretation of 
the physical behaviour of the model. 



II. THE LONG RANGE P-SPIN MODEL 

The long range p-spin glass model is defined by the following Hamiltonian 
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npi{s}) 



{l<ii<i2<---<ip<N) 



where h is an external magnetic field and the interactions Jii,i2 
distributed with the following Gaussian distribution 



P{JiiA2,---.ip) 
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exp 



are random variables 

\'-'ll,l2:---,^p I 



(2) 



The scaling of the variance with NP'^^ ensures the free energy to be extensive. For 
p — 2 the interactions are the more familiar Jij interactions and the model is the well known 
SK model For p > 2 the model is generally addressed to as the 'p-spin model'. The 
interactions between two spins Sk and si depend on all the Ji^^i2,---,ip having k and / as one 
of the arguments as well as all the remaining spin variables of the system. 

The behaviour of the model for p = 2 and for p > 2 is quite different at low temperatures. 
In this work we shall focus on the physics of the p-spin for p > 2. Note that (|l|) describes 
a model which is intrinsically mean field for large N since all the N spin variables interact 
reciprocally and there is no geometry in space. 

To solve the model, one can introduce n replicas of the system and calculate the repli- 
cated partition function, Q 



exp/3E 



E •^i.-..4---< + /^E< 



2l<---<»p 



(3) 



From now on we shall always set J = 1. Integrating over the disorder, and introducing 
n{n — 1) auxiliary fields Qab and Aafc (the Qab and \ab are symmetric n x n matrices with 
zero diagonal) one obtains the following expression 



n ^Q-" / n [-NG{Qab, Xab 



with 



G{Qab, Aafc) 
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T.Qab+lT.^-bQab-lnZlX], 



(4) 



(5) 



a^b 



a^b 



where 



Z[X] = Tr{s^} exp 



a^b 



(6) 



In equation (|6|) the spin variables have lost their dependency from the index i,-. The 
fields Xab have been introduced as Lagrange multipliers to impose the conditions 



Qab 



1 ^ 

]^E<' 



(7) 
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The saddle point of the functional integral (|4|) gives the mean field equations for the 
fields. Condition ^ gives physical meaning to the solution Qab that extremizes G{Qab, ^ab)- 
The saddle point equations for the n{n— 1) fields are 

'^ab — : 



Qab — < SaSb > — < SaSb >H{Q,S)j (8) 

where < • >h{q,S) stands for the average taken with the measure exp[—(3H{Q,S)]. 

If we indicate with and A^^ the solutions of equations (||) the saddle point free 
energy of the system is given by the following formula 



FW = + lim ^G(g,p, Xsp). (9) 



A. Replica Symmetric Solution 



In the high temperature phase one chooses the RS ansatz for the saddle point matrix, 
so one has 



Qab — q, 

Xab = A 
Qaa — Aaa 



for a b, 
for a b, 

: 0. 



(10) 



The saddle point equations deriving from this ansatz are 

q = J Dzta,nh'^{z^/X + l3h), 
where we used the following notation 



Dz = 



V2 



2tt 



■dz. 



(11) 
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In the high temperature phase one has that q = if h — 0. We shall always consider 
the h = case for simplicity of notations and because we shall not be concerned with the 
effects of a magnetic field. 

In this case one finds 



Frs 



-N 



In 2 



(13) 



This value of the free energy coincides with the 'annealed' result i.e. what would be obtained 
by calculating In Z. From (nSf) one can derive the following expression for the entropy 



Srs = N 



In 2 
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(14) 
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The entropy in (tlA) becomes negative for T < \/2\/\n2 T*, so at a greater or 
equal temperature Tc{p) the replica symmetric solution will stop to be correct. There will 
be another solution dominating Z" in the zero-n limit, and this will involve the breaking 
of the original symmetry among the replicas. The RS solution though, will not become 
unstable (if p > 2) at Tc{p) so the new solution shall be distant from the old one and the 
transition will therefore be discontinuous in the order parameter. 



B. IRSB Solution 

Here we show the IRSB solution for the p-spin model for general p. For p > 2 this 
solution is correct below Tc{p). For finite p there is a second critical temperature T*{p) at 
which the system undergoes to a continuous full replica symmetry breaking. Here we will 
not study this second transition. So one has, for p> 2, 

< T*ip) < T,{p) < 1. 

In the large-p limit one has that r*(oo) and Tc{oo) -> l/(2Vln2). The IRSB 
solution gives the following expression for Ggp 

Mm ^ = ^ {mq^ + (1 - m)ql) - ^ (mgoAo + (1 - m)qiXi) + ^ 

-i- J D{z)\n I D{y){2cosh[^„z + ^X^^„y]r, (15) 
and the saddle point equations 

A. 

qi = Dz Dy tanh^ (z + yy^Xi - Xo) cosh"^ {z + y \i - Aq) 



^ ^{q'^-q^o)i^~P) = ^ J Dz\n J Dy cosh"' (z^o + y^X^^)- 

J Dycosh"'{z^/>^ + yy/Xl - Ap) lncosh(z^/A^ + ;/VAi - Ap) 

J DycosirizV\i + yVXi-Xo) ' ^ ' 

In absence of external magnetic field qp = 0. m is an additional parameter which 
characterizes the form of the solution. The saddle point solution msp{T) has a physical 
meaning in the replica method. Two different replicas have an overlap go with probability 
rrispiT) and overlap qi with probability l — msp{T). For < msp{T) < 1 the IRSB solution 
is the correct one and the critical temperature is given by the equation 

msp{Tc) = 1. 

In the SK model {p = 2) the RS solution becomes unstable at a higher temperature 
(j^SK _ 2 with our normalizations) and there is a transition into a full RSB phase with a 
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continuous solution q{x) which holds for all temperatures below T^^ = 1. In this case the 
IRSB is anyhow a better approximation than the RS one. 

In figure (|l|) we plotted the free energy of the model for p = 4 for low temperatures. 
The high and low temperature curves are tangent at the critical temperature. In figures 
we plotted the values of qi and m obtained by the IRSB solution for p = 4 in function 
of the temperature. 
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FIG. 1. Free energy of the p-spin model vs temperature for p 
solution which is wrong at low temperatures. 
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The dashed line is the RS 



0.93 I 

0.25 



0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 

FIG. 2. qi vs temperature for p = 4. 
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FIG. 3. m vs temperature for p = 4. At the critical temperature the curve crosses the line m = 1. 



III. THE SHORT-RANGE MODEL 



Here we introduce and study a short-range version of the p-spin glass. The model that 
we introduce is defined on an hyper-cubic d-dimensional lattice of side L. On each side of 
the lattice there arc M spins. Every spin interacts via quenched random couplings with 1 
spins chosen among spins on the same site and on nearest-neighbour sites. It is natural to 
expect (and it will also derive from the equations) that, for large M, one recovers the mean 
field solution since each spin interacts with a large number of other spins. The Hamiltonian 
of the model is 



M 



(17) 



0^ ii )" ■ 



By X]<ji I > snuY over all the sites of the lattice taking, for each couple; of adjacent 
sites i and j, p — l: of the h, - ■ ■ ,lp indices equal to i and k indices equal to j. In other words, 
for each nearest neighbour sites i and j, every interaction involves p — k spins of site i and k 
spins of site j with k running from zero to p. We consider discrete (±1) spin variables and 

we call .Si^, (Z, ) the i*'' spin of site Ir with ir running from 1 to M. The j]^ are quenched 
random variables extracted from the distribution 



aMP- 



■ exp 



(18) 



where the normalization ensures an extensive free energy and the constant a will be fixed 
later on by the condition that, in the large-M limit, the free energy density of this model 
coincides with the one of the long range model. 

The model just defined is a short range model, but the existence of the parameter M 
ensures the mean field solution to be correct in the limit of M ^ oo. In this way one has 
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a systematic way to reduce arbitrarily the corrections to mean field which are due to the 
finiteness of the coordination number and one can control the crossover effects from finite- 
dimensional to mean field behaviour. A similar generalization of the Random Energy Model 
has been studied in the one dimensional case in | |lO|] . 

The mean field solution of the model is identical to the solution of the long range p-spin 
model. Nevertheless, even in the large-M limit, this model keeps an important ingredient of 
difference with respect to the long-range p-spin model because there is a geometry in space. 

The expression for the partition function of n replicas of the system is 



a<b I 



n n / dQab{l)dXabil) exp [-MG[Q, A]] 



(19) 



with 



G[Q,A] 



'2 " 



2a 



In Tr exp 



a<b I 



E E {Q'ab---Q'ab) +T.T.^abQL 

a<b <li,- -,lp> 

n L"^ 

EE^^^-^^w^'w 

a<b I 



(20) 



We can write the non local part of the G[Q, A] in a more convenient form exploiting the 
following identity 



(21) 



l.m 



where 



and ^he sum over all possible directions of the unitary lattice spacing vector 1. So 

the function G[Q, A] takes the form 

G[Q, A] - E E ^('' (Q'-b + QabY - (2^ + 2d - 1) j2{QUr 



a<b l,m 
L" 



EE^ahQih-lnTrcxp 



a<b I 



a<b I 



(22) 



The mean field solution must be constant in space because the first term in ( po|) can be 
re- written as a gradient plus a local term. We want to normalize the G[(5, A] in such a way 
to recover (^ in the mean field approximation. To do so one has to set 

a = 2Pd-2d+ 1. 



The saddle point equations yield 
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Qi,^{S'^il)S''il))^. (23) 

These equations are the same equations that we found in the long range case and are 
exact in the hmit of M — > oo. Therefore, in the large-A/ hmit, the equations (|23| ) reduce to 
equations ( [n] ) or jl^ ) depending on the temperature. 



IV. GAUSSIAN FLUCTUATIONS 



For finite M the saddle point equations are not exact. In this section we estimate the 
0{1/M) corrections to the free energy due to the small fluctuations around the saddle point 
of both the high and low temperature phase. For finite M we can write 



Qab ^ Qal + ^Qab ^ab = + JA^^, 



(24) 



where 5Q^^^ and 5\'■^^^ are 0{1/M) quantities. and A^^ are the solutions of the equation 
(23) and are site independent. 

It is convenient to work in Fourier space so we write 



SQL 



. (2^)' 



<5Qa6(fc)e'" 



k-l 



(25) 



One can expand G[Q, A] up to quadratic order in 5Q\^ and (5A^j,. Observing that the 
couple (a/3) takes n{n — l)/2 different values, we can write 



2 J-^ ^^-y 



naik)Maf3np{-k) 



(26) 



where MFsp is the free energy of the system in the mean field approximation. We have used 
the following labeling: 



ria(fc) = SXab{k) for 
ria(fc) = 6Qab{k) for 



1 < a < 
n(n — 1) 



n{n — 1) 



< a < n{n — 1). 



Map is the n{n — 1) x n{n — 1) matrix of the fluctuations 

M = 



"^abcd ^^^abcd 
M<}'^ Afl'l 
^"abcd ^"abcd 



(27) 



(28) 



with 



'■^'-abcd 



/3Mp-1) ^.p(p-2), l-2d + 2Pcos^(|) 
2 ^'^'^ "'^"'^ \-2d+2Pd 



abed [{S^S'^S^S'')^ {S^S')^ {S^S'')^ 



Mabcd 



bed- 



(29) 
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We indicated with Cabcd the connected four-spin correlation function. We also used the 
notation 



Performing the Gaussian integral one obtains 



F{T) = Fsp{T) + ^^FiT), (30) 



with 



= ^ 1™ E r 1^ [^(k)] , (31) 

where m\ is the multiplicity of the X{k) eigenvalue of Map. 

If the temperature is above the critical temperature the mean field solution is RS with 
qo — Q (we shall always consider the case of zero magnetic field). The matrix of the fluctu- 
ations has got the whole block M^l^^ equal to zero and it is straightforward to show that 
above the temperature-dependent part of the 0{1/M) corrections is zero and 

AF = 0. (32) 

In the low temperature phase the structure of M^p changes because the sub-matrix 
■^abcd ^'-'^ identical to zero but takes the following form 

^abcd = if a and b do not belong to the same block 
^abcd ~ fi^^ ^) ^ ^ belong to the same block , 

where 

/(,, k) = _PMP-l)^,-.l-^d+^^cos^^,) 

■'^^^ ' 2 ^ l-2d + 2Pd ^ ' 

In this case one needs to calculate all the n(n — 1) eigenvalues of Ma^p- The calculation 
is similar to the one performed in |l^ in the case of the Little model with the difference that 
in that case the model was Gaussian and there was no need to introduce the Xap field. So 
here the eigenvalues are n(n — 1). We shall indicate with Aa and Ag the eigenvalues of the 
two sub-blocks M^^^^ and Af^^^^. It is easy to calculate the determinant of the total matrix 
M observing that, because of the diagonal structure of A/^^^^, correspondingly to every Aa 
and Xq there are two eigenvalues of the total matrix 



Xa + Xq ± ^ (Aa — A, 



l2 



A± = - — " — — ' 

so that A+A_ = AaAq — 1. 

The structure of the spectrum of the eigenvalues Aa and Xq in the IRSB phase is briefly 
reported in the appendix. The final expression for the O(-jy) corrections below becomes 
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Im 

m 

m — 1 
m 



1 r d'^k r(m-3) 



2/3 {2TT)d 



■ln[l + /(q,k) {C[abcd] + il-qf)] 



i-/(g,k) 



A + D + \A-D\ 



In 



l-/(9,k) 



A' + D' + M' - D' 



(35) 



where we have defined 



/ cia rib ric / C^I /c^Cf^X 

L^[a6cd] — W -J /a ^ W '-' /A W /A ' 

with the four repHcas a, 6, c, befongingto the same diagonal block of the IRSB matrix. 
The constants A, D, A' , D' in equation ( p5[) are also defined in the appendix. 




-0.838 

0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 

FIG. 4. Free energy density of the three-dimensional p-spin model vs temperature for p = 4. and 
M = 4,6,8, 10, oo 



As a check one can easily verify that, for m = 1, the result (32) is, as it should be, 
recovered since the structure of the equations becomes identical to the RS case after the 
identification qi = qo- 

To perform the above integrals in momentum space we make use of the identity 



d'^k 
.(2^ 



In 



l + Acos2(-) 



/ — cxphi] 
Jo ^ 



At \ ^ ( Adt\ 



(36) 



where /o is the modified Bessel function of order 0. 

One can solve numerically the saddle point equations and then plug the value of 
qi{T), m{T) into the expression of G and AG. The result, for d — i and M = 4, 6, 8, 10, oo, 
is shown in figure (^). In the plot the higher are the curves the smaller is the value of M. 
We can see that the perturbative corrections slightly shift the curves but do not change the 
nature of the transition. The critical temperature, defined as the temperature where the 
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curves of the high and low temperature free energy coincide, is the same at which rUsp = 1. 
Furthermore, the two curves are tangent at the critical temperature denoting a zero latent 
heat as in the MF case. Finally, we also remark that the transition temperature changes 
very slightly with M. A rough estimate the 0{1/M) shift in the critical temperature gives 

T,(M)~T,- ^(0.0168). 



V. REPLICA PROPAGATORS 
A. At equilibrium 

The study of the Gaussian propagators provides knowledge on the existence of zero- 
mass modes and therefore on the nature of the transition. In this section we consider the 
corrections to the mean field solution of the theory defined by the (p2|). In a perturbative 
approach we define and calculate the propagators of the p-spin model in the RS and IRSB 
phase. 

The propagators Gated are defined in the following way 

Gabcd- (SQabSQcd) 

Gtbcd^ (SAabSA^d). 
The equations for the propagators are easily found by standard arguments 

and in the case of our functional integral ( pj| ) they become 

dG[Q,A] ^ dG[Q,A] ^ 

{oQab—^ ) = Oab.cd \"^ab—^ ) = 

oQcd oQcd 

dG[Q,A] ^ dG[Q,A] ^ 

{oQab TTl ) = {oAab TTl ) = Oab.cd- (38) 

oAcd oAcd 

It is straightforward to obtain the following linear equations for the propagators 
2f{Q,d,k)Gt^^^ + 2G^dab^0 

"^Gabcd ~ 2 ^ CrscdG^^ib = 
r<s 

■^Gabcd ~ 2 5Z GrscdG^sab — ^ab,cd, (39) 
r<s 

where /(Qcd, k) is defined by (|3^). 

We shall always use the fact that, in absence of magnetic field, qo — 0. There are 
different kind of propagators depending on the relation between the replica indices they 
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refer to, so for different choices of indices abed in equations 
equations. In the RS phase one obtains 



^abcd — 



'-'abed ~ 



2 ^ab,cd 



'-'abrd ~ n^abcd, 



^) we will obtain different 

(40) 
(41) 

(42) 



where Cabcd is the four-spin connected correlation function defined in (p9|). So in the replica- 
symmetric phase all propagators G^j,^^ are zero except for the diagonal one which is trivial 

(43) 



G 



bab 



G 



hba 



In the IRSB phase, we shall use the convention to indicate with [ab] two replicas be- 
longing to the same block and with [a\\b] two replicas belonging to different blocks of the 
IRSB saddle point matrix Q*^. If a and b (or c and d) do not belong to the same block , 
the form of the propagators is still given by equation ( 42 ) where the function Cabcd changes 
in the low temperature phase so one has 

1 

2' 

1 r , 





- rQ 


[b] 


a][b] 


rQ 


- rQ 


[b] 


c][b] 


rQ 


- rQ 


[b] 


c][d] 



91, 



M. 



(44) 



where a and b belong respectively to the same blocks as c and d. 

For the replica indices in the same block, [a6], we have Qab = qi and from equations 
(39) we obtain the following equations for G^ and G*^ 



= /(<zi, k)-i + fiqu^r'Gt 



Gabcd - -/(^l, k) CrscdG. 



A 

/ J ~ ' - rsab 



bed 

/(gi,k) 



b.cd- 



(45) 
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FIG. 5. Gi , Ga , G* for p = 4 and d = 3 

The second of the equations (^) can be written for the various choices of the indices 
abed obtaining a set of four coupled equations for Gi, G2, G3, G4 defined as fohows: 

Gl ~ G[af,afc] G2 — G\ahar\ 

Gz = Gyahrs\ Gi = G[ab][rs] (46) 

Solving the equations and plugging in the MF values for qi and m we obtain the 
propagators in function of the temperature. From the first of the equations ( |45| ) it is possible 
to define the propagators G^j,^^. One has that, in absence of magnetic field G4 (fc = 0) = 0. 
In figure (||) we plot the Gi,G2,G^ for fc = versus the temperature. In conclusion, we 
observe no divergence of the RS and IRSB propagators at the critical temperature. 

This means that no zero-mass modes are present around the stable saddle point solutions 
and usually it indicates that no continuous transition is taking place. 

B. On the dynamical line 

We recall that in the long range p-spin model, there is a dynamical critical temperature 
Til > Tc below which the system does not reach equilibrium in finite times. In the mean field 
approach, if we quench the system to low temperature coming from an high temperature 
state, the system goes to a dynamical metastable state, having an energy greater than the 
equilibrium one. In |]l3|] , it has been shown that if the system starts from a random initial 
configuration, it evolves following the flat directions connecting an ensemble of metastable 
states which are the threshold solutions to the TAP equations jlj]. 

This cannot happen in short range models and it is and artifact of mean field theory. 
However it can be interpreted as the signal of starting a quite slow approach to equilibrium. 
In any case the disappearance of this metastable energy is a non perturbative result, so that 
it make sense to study the properties of the correlation functions in this state. 
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In the replica approach it can be shown that in the dynamical metastable state the 
value of m is not fixed my extremizing the free energy, but by imposing the condition that 
the replicon eigenvalue is equal to zero. We are going to compute the propagators in this 
region as functions of the temperature. The formulae are similar to the previous ones, only 
the value of m will be different. 

This amounts to set 

/(g,k = 0) = ^ (47) 

where is defined in the appendix. 

We thus obtain some new equations for the propagators concerning replicas belonging 
to the same block 

_ (6 - 5to + TO^)(g - ro)(4g - mg^ - Arg + mrp) 

A{f{q,-k)-lo){f{q,-k)-lo-{m-l)h)if{q,k)-h) ^ ' 

^ (3 ~m){q- ro)(4g - mq^ ~ 4ro + mrp) 

A(/(gk,)-M(/(q,k)-/o-(m-l)/i)(/(gk,)-/2) ^ ^ 

^ ^ 2(g - r-o)(4g - mq^ - 4ro + mrp) 

' A(/((z,k)-lo)(/(9,k)-Zo-(m-l)Zi)(/(gk,)-;2) ^ ' 



where 



A = (-1 + 2(7 - ro)^(-2 + IQq - Qmq - i2q^ + 2'im^q - im'^q'^ + Amq^ - dm'^q^ + 
m^q^ — 12ro + 7mr0 — ni^ro + 4:8qro — Admqro + ISm^qro + 
—m'^qrO — Zmq^VQ + 4m^g^ro — w?q^rQ — ISrg + 21mr^ — Srri^rQ + m'^rg) 
-1 



l-2q + ro 

1 



^0 

h 



-1 + 4(7 — m(7 — 3ro + mro 
-2 

2 — 8(7 + Amq + 771(7^ — 777.2(72 + 6ro — Smro + 7772ro 



In the previous formulae, rg — C^abcd] ■ In the general case the propagators diverge, for 
small k, as k^^. The condition m = 1 now gives the dynamical critical temperature, and 
one can see that all the propagators in ( pO| ) coincide and that the divergence is of order k~^. 

The divergence at small momenta of the propagator is not unexpected . Indeed the 
vanishing of the replicon propagator implies a divergence of the susceptibility and conse- 
quently a singularity at = 0. The form of the singularity could not be predicted using 
general arguments. The change of the exponent at Td is particularly striking; it is related 
to the degeneracy of the replicon and longitudinal eigenvalue at to = 1. 

Although the result was derived in the context of spin models, we believe that this 
structure of exponents fc^^ at Td and fc~^ at T < is quite general and it would be valid in 
many others models. It is also clear that these are mean field results: also if we remain in a 
perturbative framework, the exponents are likely to be changed for sufficient low dimension. 
The value of the upper critical dimension (6?, 8?), above which the mean fields exponents do 
not get perturbative corrections, can be extracted by analyzing the contribution of higher 
loops, but this task goes beyond the aim of this paper. 



15 



VI. DISCUSSION 



The model that we have described has been studied numericaUy in d = 3 with p = 4 and 
M — 3,A for T > Tc and below Tc The results of the numerical simulations seem to 
indicate, the existence of a transition at a critical temperature Tc from a high temperature 
phase to a broken replica symmetry phase. However, this transition appears to be of second 
order with divergent spin-glass susceptibility. A possible interpretation of this apparent 



contradictory phenomenology is the following |17|. For each realization of the disorder there 



are some regions of space in which the effect of frustrations are weaker than in the rest of 
the system. Within these regions the system is likely to freeze at a temperature Tr which is 
higher than the temperature Tc at which the whole system freezes. So for T > Tc (not too 
high temperatures though!) there are regions in space which the system is locally frozen, 
the typical size of the region being a function of the temperature i?'^'' (T) where dr is the 
dimensionality of the region. Within these regions the system is very strongly correlated 
and the correlation length is of order R{T). The total SG susceptibility is the integral over 
space of the local SG susceptibilities and the contribution of the regions in which the system 
is strongly correlated grows with R{T) and diverges when R{T) becomes of the size of the 
whole system. So, in this interpretation, the transition remains discontinuous within the 
regions of space where it occurs, i.e. the local overlap changes discontinuously from qo to 
qi , and the continuous varying quantity is the typical size of the regions where the system 
is strongly correlated. This work wants to be a step forward towards the comprehension of 
spin glass models and of structural glasses in finite dimensions parallel to other attempts 
There are still many unclear things on the subject which is worth, according to us, 
for further studies. 
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APPENDIX A: 

In the IRSB phase there are three invariant subspaces of eigenvectors corresponding to 
three classes of eigenvalues. In this appendix we shall list the eigenvalues \\ and Aq and 
their multiplicity fi for each different subspace. 



1. Longitudinal eigenvalues 

For each diagonal sub-matrix there are two couples of longitudinal eigenvalues with 
multiplicity jiL — ^ 

a;- = ^±^4^ = /(,., k, 

where 
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„2 j_ _ „2^ (to - 2)(to - 3) ^ 

2 

£>=-l + 2(m-l)gi + (m-l)^g^ (A2) 



A= 1 - g2 + 2(m - 2)(<zi - g?) - V'" (Al) 



2. Anomalous eigenvalues 

For each diagonal sub-matrix there are four couples of anomalous eigenvalues. 
Two of them have multiplicity /z^ = (n — m)/m and coincide with the longitudinal 
ones. 

The other two have multiplicity = n{m — l)/m and are 
where 

A'= l-ql + {m^ 4)(5i - ql) - (m - ^)C^abcd] (A3) 
D'= -l + {m-2)qi- {m-l)ql (A4) 

3. Replicon eigenvalues 

There are four couples of different replicon eigenvalues 



= Pi - 2Qi + i?o A0 = /(gi,k), = n^^^^ 

A^ = Po + 2(m-l)Oo + (m-l)V = 0, Mfl = n ^"~Jt"'^ 

Af = Po + (m-2)Qo-(m-l)g? Ag = 0, = n ^" ~ ^"^^^"^ ~ 



(A5) 
(A6) 



where 



Q 

(A7) 

Af = Po-2Qo + 9? Ag = 0, Mi^ = n ^"~";^^7~^^' (A8) 



Po = -1 Pi=ql-l 

(A9) 

Qo = -91 Ql = 9i - 91 

(AlO) 

Ro = Ql - C[abcd] (All) 
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